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We study the potential to use measurements of the properties of extra neutral gauge bosons 
(Z''s) in pp collisions at the Large Hadron Collider to unravel the underlying physics. We focus 
on the usefulness of third generation final states (r, h, t) in distinguishing between models with 
non- universal Z'-fermion couplings. We present an update of discovery limits of Z''s including 
the 2010-2011 LHC run and include models with non-universal couplings. We show how ratios of 
o{pp — 5> — > ti), a{pp Z' ^ bb), and a{pp — 5> — ^ r+r") to a(pp Z' ^ M^^^M") can be 
used to distinguish between models and measure parameters of the models. Of specific interest are 
models with preferential couplings such as models with generation dependent couplings. We also 
find that forward-backward asymmetry measurements with third generation fermions in the final 
state could provide important input to understanding the nature of the Z' . Understanding detector 
resolution and efficiencies will be crucial for extracting results. 

PACS numbers: 12.60.Cn, 12.15Ji, 14.70. Pw, 12.15.-y 



I. INTRODUCTION 

With the startup of the Large Hadron Collider (LHC), 
it is possible that direct experimental evidence of new 
physics will soon follow. Many theoretical models of 
physics beyond the Standard Model include a new colour- 
less, spin-1, neutral gauge boson (Z') [J-Q. A kinemat- 
ically accessible Z' would be one of the most distinctive 
signals at the LHC, and may potentially be one of the 
earliest discoveries. If a heavy neutral gauge boson were 
discovered, the immediate task would be to understand 
its origins by measuring its properties. This subject has 
been explored extensively and continues to be an impor- 
tant topic of on-going research 6-14]. 

The ATLAS pTs] and CMS [iGj experiments at the LHC 
expect to be able to identify r leptons and h- and t-quark 
jets in pp collisions. In a recent letter, we suggested that 
third generation fermions can provide an important tool 
for discriminating between candidate models of Z"s 
Along this vein, we found that such measurements can 
also be used to explore the gauge symmetry of some 
models of new physics that include a Z' with distinc- 
tive properties reflecting the mixing between the under- 
lying gauge groups. These include models with fermion 
type dependent couplings (i.e., preference for quarks over 
leptons) p7l - IT9| . or generation dependent couplings (i.e., 
preference for third generation) [2o!-[3l|. The models 
with generation dependent couplings, such as extended 
technicolour [24l - l26| and topcolor assisted technicolour 
[l^-dll) are especially interesting as they take the heavy 
top quark mass as evidence for new physics that is con- 
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nected to the mechanism of electroweak symmetry break- 
ing (EWSB). Thus, evidence for the violation of genera- 
tion universality could give important hints to the mecha- 
nism of EWSB. For such models, standard analysis meth- 
ods that rely on Drell-Yan processes with either electron 
or muon final states may fail to identify important details 
of the underlying physics. 

In this paper, we look at a number of measurements 
involving third generation final states, expanding on our 
previous letter [7]. We start by exploring the discovery 
reach for the models we consider. However, our main fo- 
cus is to consider the viability and usefulness of measure- 
ments involving third generation fermions in identifying 
the underlying theory of a newly discovered Z' . The 
measurements we study are the ratio of the third gener- 
ation quark cross sections to the muon cross section, the 
ratio of the tau cross section to the muon cross section, 
and the forward-backward asymmetry with bottom and 
top-quark final states. 

We found that the ratios of the branching fraction of 
T^r~, hh and tt final states to the [i^ [i~ final state in 
Z' decays are especially useful for understanding models 
with preferential couplings and can be used to quantify 
the extent of the preference. We also found that ApB 
measurements with third generation quarks in the final 
state have the potential to provide Z'-fermion coupling 
information. Several papers have been published in the 
past on related topics |32| - |4q |. This paper takes into ac- 
count recent information on the expected capability of 
the ATLAS [15,] and CMS J^] detectors in addition to 
backgrounds. 

We begin in Section II with a brief survey of models 
with non-universal couplings to ^"s, followed in Section 
III with details of our calculations, focusing on issues 
of third generation fermion identification and Standard 
Model backgrounds. In Section IV we give our results. 
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starting with an update of discovery limits for Z"s from 
the models we consider in this paper and reflecting the 
updated LHC running plans. Our main focus, however, 
is exploring how one can use third generation fermions to 
learn about the underlying theory that gives rise to a Z'. 
Finally, we summarize our main conclusions in Section 
V. 



II. MODELS 

New massive resonances are present in many models, 
including KK theories with finite size extra dimensions, 
string theories and theories with extended gauge sectors. 
In this paper we focus on models with extended gauge 
sectors. Some examples of these models are variations of 
Little Higgs Models (Littlest Higgs (LH) with tan (j)H — 
1.0 [ill. Simplest Little Higgs (SLH) [42|, Anomaly Free 
Simple Little Higgs (AFSLH) @), Eq models (x, v), 
Left Right Symmetric models (LRM, ALRM) {ga = .g^), 
and 3-3-1 models [Hi. The details of the E& and LR 
Symmetric models have been described elsewhere, so we 
refer the interested reader to the literature 

We single out a subset of models where spontaneous 
breaking to diagonal subgroups results in generators cor- 
responding to gauge bosons that couple to either differ- 
ent generations (first vs third) or different types (quark 
vs lepton) of fermions with differing strengths. These 
models typically have a group structure of the form 
SU{2) X SU{2) or U{1) x f/(l). Such models include 
the Ununified Standard Model [l7|, no n-commuting ex- 
tended technicolor [23] , top-flavour I2l| - [23| , and topcolor 
assisted technicolor models [27]. These models can be 
distinguished using third generation fermions and we 
therefore give some relevant details of the models that 
are used in our analysis. 



general, Mz' > 2 TeV is required, with a stronger mass 
constraint for larger values of sin^ 4>uum- In our study, 
the use of Mz' = 1.5 TeV is for comparison to other mod- 
els without such constraints. Unless otherwise stated, we 
take sin, (tjjnj m =0.5, following the work of Chivukula et 
al. 

B. SU(2)h X SU(2)i - Extended Technicolor (ETC) 

In models of an extended SU{2)h x SU{2)i struc- 
ture, fermion generations transform differently under 
each gauge group - the first two generations transform as 
a doublet under the SU{2)i, while the third generation 
transforms as a doublet under the SU{2)h- Subsequently, 
this extended gauge group is broken to its familiar diag- 
onal subgroup, SU{2)l, at some energy scale fi. The 
electric charge operator is given by: 

Q = ni + n,, + Y. (2) 

In diagonalizing the mass matrix for the neutral gauge 
bosons, a nearly SM Z° arises with a Z' that couples 
only to left handed fermions. This Z' then has a greatly 
enhanced coupling strength to third generation fermions, 
with couplings given by: 

9l = gzocyj (- tan (jjETcTsh + cot c/jetcTsi) . (3) 

Unless otherwise stated, a value of sin^^TC = 0.9 is 
used. Constraints found by Chivukula et ai, [l^, [25| 
require Mz' > 2 TeV for this value of sm(j>ETC- Again, 
calculations are performed with lower values of the Mz' 
for comparisons with other models. 

Examples of these models include top-flavour mod- 
els [2l| - |23j and non-commuting extended technicolor [23 - 



A. Ununified Model (UUM) 

In the Ununifled Model the left handed quarks 
and leptons transform as doublets under their respective 
SU{2) groups in SU{2)q x SU{2)i x U{1)y. Right handed 
fermions transform as singlets under both groups, and 
hypercharge assignments remain the same as in the SM. 

After symmetry breaking, the mass eigenstates for the 
gauge bosons include a massless photon, and two mas- 
sive Z bosons. The light Z° boson deviates from SM 
couplings at order sin^ cj), where (j) is the mixing parame- 
ter between the two SU{2) groups and is expected to be 
small. The heavy boson then couples to the lepton and 
quark sectors as: 

gz' = gz«Cw f 7 — ta.n (jyuuMTsi) . (1) 

V tan 0(7(7 A/ J 

Chivukula et al. [H, [l^ found 95% C.L. constraints on 
the value of Mz' dependent on the mixing angle. In 



C. U{l)h X I7(l)i - Topcolor Assisted Technicolor 
(TC2) 

The large mass of the top quark has led to suggestions 
that the top quark is intrinsically related to the dynam- 
ics of electroweak symmetry breaking. Models based on 
this idea are often called Topcolor [4^ , and contain 
a gauge structure that is generation dependent. It is as- 
sumed that the QCD gauge group and the hypercharge 
C/(l) arises from the breaking of a larger group as in 
S'C/(3)i X SU{3)2 X U{l)i X U{1)2 ^ SUi3)c x U{1)y 
with a residual, high scale S'[/(3)' x J7(l)'. In this type 
of model, the Z' plays a role in the generation of the large 
top quark mass by providing a tilting mechanism for the 
top quark seesaw [35l |. 

The third generation fermions transform under the 
C/(l)i group, while the first two generation fermions 
transform under the U{1)2 group. After symmetry break- 
ing, the resulting Z' couples differently to the third gen- 
eration fermions than it does to the first two generation 
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fermions. The couplings of the SM fermions to the Z' 
are given by: 

3 1 

9l,R = -j9Z°S^YsM cot (/)TC2 (4) 
1 2 1 

9l,r= 29z<'SwYsM'^'d''a.<i)TCi-, (5) 

where Ysm is the standard model hypercharge value. Un- 
less otherwise stated, a value of sok^tci = 0.5 is used. 



and 



In Eq. ([6]) through ([9]), Mi and Vi are the masses 
and widths of the photon, SM Z°, and Z'\ L^, i?}. 
are the left and right handed couplings of the gauge 
bosons to fermion species /; rrij is the mass of the fi- 
nal state fermion; fq^q{x,Q^) are the parton distribution 
functions (pdf's); Xa and Xb are the momentum frac- 
tions of the partons; is the scale at which the par- 
ton distribution functions are evaluated, which we take 
to be s - the square of the parton centre-of-mass en- 
ergy; 9 is the centre-of-mass scattering angle; and (3/ = 

^1 — Airij/s. In our calculations we used a = 1/127.9, 

sin^e^ = 0.231, Mz = 91.188 GeV, Tz = 2.495 GeV 
and nit = 171.2 GeV [i^, and for the pdf's we used 
set CTEQ6M We calculated the Z' width includ- 

ing only decays to SM fer mions, and neglected 

Z — Z' mixing and decays to pairs of W and Z° bosons, 
as well as decays to heavy, fourth generation fermions. 

The total cross section is proportional to the Sq and 
S'^ terms (with summed couplings), which are depen- 
dent on symmetric combinations of the Z'-fermion cou- 
plings. The antisymmetric combinations of the Z' - 
fermion couplings in the Aq term contributes to the 
forward-backward asymmetry, ApB- Thus, measur- 



Examples of these models also include top quar k see- 
saw models [27n29iJ and flavour- universal TC2 |30l, [3l|. 



III. CALCULATIONS 

For our calculations we use the leading order Drell-Yan 
cross section. While it can be found in the literature, we 
include it here for completeness [1, [13, Ell : 



ing the production cross section and forward-backward 
asymmetry of the decay of a Z' to third generation 
fermions will give complementary information about the 
Z' couplings to fermions - information that will be cru- 
cial for disentangling the underlying theory. The mass 
dependent coefficient for 5*^ is only relevant for top de- 
cays when the mass of the Z' is relatively small. The 
effect is 0(10%) for a 1 TeV Z', dropping to 0(1%) for 
a 2 TeV Z\ for top quark final states. 

In our calculations, we include a K-factor [51] in the 
cross section to account for NLO QCD corrections and 
neglect NNLO as the uncertainties in parton distribu- 
tions dominate over such small effects [H, [l^ . We also 
included QCD [Hi-IHi] and EW (ES-Hi corrections in the 
width, with weak corrections having only a very minor 
effect. 

The formula given in Eq. |6] was used to calcu- 
late the cross sections and distributions in this pa- 
per. The phase space integrals were performed using 
Monte-Carlo integration methods with weighted events. 
QCD backgrounds are calculated using a combination 
of the WHiZaRD (with O'MEGA matrix elements) 
(sl and MADGRAPH [63] Monte-Carlo event generators 
and compared with our own code following Barger and 



'^'^ ^pp^Z',Z",7^//)-V / rfx,dx,/,(x„Q^)/,(:.,,Q^)^ + /,(x„Q^)/,(x,,g^) ^'^j^ /\ (6) 

„ J a cost' dcosO 



d cos 6 ^ 



;re da /d cos 9 is given by 



dcos( 



qq^Z\Z\-l^ ff) = ( (l + /3/ cos^ 9) Sq + 2/3 f cos 9Aq + S'A (7) 



r 
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Phillips [63]. All three background calculations agreed. 



IV. IDENTIFICATION AND BACKGROUNDS 

The challenge to using third generation fermions in 
LHC events is to first identify the third generation 
fermions and then to identify the signal events buried 
in the large Standard Model QCD backgrounds. The 
backgrounds take two forms. The first is the SM pro- 
duction of third generation quarks which we will discuss 
below. The second is the misidentification of the large 
QCD light jet backgrounds as heavy fermions. Identifi- 
cation of third generation fermions and misidentification 
of light jet backgrounds are not unrelated. The rejec- 
tion of light jet backgrounds is crucial to making accu- 
rate measurements of heavy quark final states. However, 
higher rejection typically results in lower identification ef- 
ficiencies, and correspondingly lower statistics. There is 
therefore a tradeoff between high identification efficiency 
and suppression of the reducible background. 

There has been recent activity on this subject with a 
number of papers appearing on the topic of the identifi- 
cation of top jets, specifically with techniques for highly 
boosted jets |3S, i65. - ;67j . In addition, ATLAS has re- 
leased updated information on the expected capabilities 
of the detector for the tagging of leptons and bottom 
quarks [l^ . We use this information to estimate the effi- 
ciency and rejection rates needed to analyze the data. It 
should be noted that some estimates are given for tagging 
lower energy events and we do not evaluate the validity 
of extending or extrapolating to higher energies. Bet- 
ter understanding of the ATLAS and CMS detectors will 
ultimately give more reliable values. Finally, there are 
non-QCD backgrounds such as W+jets, (Wb+Wb), and 
Wbb + jets. However, studies have shown that these can 
be controlled by constraints on cluster transverse mass 
and the invariant mass of jets [36| so we will not consider 
them further. 

In the following subsections we summarize the identi- 
fication efficiencies and fake rejection rates for the third 
generation fermions that we use in our analysis. 



A. Muon 

Muons are the most distinct of all tagged signatures for 
both ATLAS and CMS |68|]. Recent studies by the AT- 
LAS Muon Working Group [l^ show a single muon iden- 
tification efficiency greater than 95% for px > 30 GeV. 
These values vary over pseudorapidity and energy; for the 
purpose of this analysis, we use an efficiency of — 96% 
for single muon identification and consider only Drell-Yan 
backgrounds. 



B. Tau Lepton 

There are three possible decay modes for a t+t~ 
event: purely leptonic {t~^t~ — )• Vr^rl^vl^v, 12.4% of 
total events), semi-leptonic events (r+r^ VtVtIv -|- 
jeis, 45.6% of total events), and purely hadronic events 
(r+T^ — > VrV-r + jcts, 42% of total events). Of the 
hadronic decays, sa 77% have one charged pion track 
(one-prong) and « 23% have three charged pion tracks 
(three-prong). In our analysis, we consider the results 
from the purely hadronic mode, which suffers from a large 
reducible dijet background, as the leptonic modes have a 
larger missing energy component. It has been shown that 
a T+T^ invariant mass distribution can be reconstructed 
in these modes even though the r decays produce missing 
energy 

In order to reduce the dijet background, it will be nec- 
essary to implement tight tagging methods. In partic- 
ular, ATLAS [15] estimates that it may be possible to 
achieve a rejection rate of 10'^ and an efficiency of 20% 
for 3 prong decays with Et > 100 GeV, and a rejection 
rate of 10^ and an efficiency of 50% for 1 prong decays 
with Et > 100 GeV. This amounts to an overall rejection 
rate of 5 x 10^ with w 20% of aU hadronic (42%) t+t~ 
events passing the selection criteria. The number of use- 
ful T'^T~ pairs can potentially be increased by including 
the leptonic and semileptonic modes. 

C. Bottom 

According to the current ATLAS algorithms [l5j] , a jet 
is tagged as a &-jet if a b- quark is found with pT > 5 GeV 
within AR = 0.3 around the centroid of the jet. Further 
improvements are made to this tag by including either the 
identification of a secondary vertex or else reconstructing 
the invariant mass of the jet. In an analysis of 6-jets 
arising from WH and tt, a jet rejection of about 2 x 10^ 
was found for a 6-jet efficiency of 60% - or 36% for a bb 
pair, with a light jet rejection of 4 x 10''. Tighter tagging 
could be employed to improve the rejection rate if the 
light jet background dominates the observed events. 

D. Top 

A complete analysis of physics events involving top 
quarks would require an analysis of the decay products to 
calculate an overall efficiency. This is beyond the scope 
of this paper and we instead use previous studies to es- 
timate an overall efficiency (36l . l65l . [7ll - l73j . 

We considered both the semi-leptonic decay mode 
(tt — > (lvi){jj)bb where / = /i, e) with an overall 
30% branching ratio, and fully hadronic modes {tt — > 
{jj){jj)bb) with an overall branching ratio of 46%. 

Semi-leptonic modes present the best signal-to- 
background ratio, with the dominant fake background 
being W + nj. An analysis of KK modes for the LHC 
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by Agashe et ai, [zl] found an overall efficiency of « 1%, 
after taking into account branching fraction, cuts and 
6-tagging efficiency. Their analysis included a 20% ef- 
ficiency for the 5-tag. They found the reducible back- 
ground to be small compared to the irreducible QCD 
background. 

In a separate analysis, Kaplan et at. [65l | considered a 
method of top tagging for high px hadronic top decays 
that suggests better efficiencies and does not incorpo- 
rate a 6-tag. Their analysis using Pythia and a detailed 
systematic examination of top decays suggests a tagging 
efficiency of approximately 35% to 45% for top jets with 
transverse momentum between 600 and 1400 GeV. For 
fully hadronic ti events, this technique could result in 
an efficiency > 10% for «46% of events, and a rejection 
against light jets > 10^. Very specifically, the dijet re- 
ducible background is shown to be reduced to the same 
level as the QCD tt background, where a resonant peak 
may be visible. 

Almeida et al. (73j found that using only a cut on the 
jet mass resulted in a single top jet identification effi- 
ciency of between 34% and 58% with a rejection factor 
of around 30. By applying further cuts on the jet struc- 
ture, they were able to achieve a rejection factor of about 
5000 for a top identification efficiency of ~2I%. This is 
in general agreement with the Kaplan study in that rea- 
sonable top jet tagging efficiencies can be expected at the 
LHC without significant light jet contamination. 

In our analysis, we assume the method used by Kaplan 
et al. and employ an efficiency of 16% (40% for each top 
jet - the midpoint of the expected range) for the fully 
hadronic decay mode, with a rejection of dijets by a factor 
of lOl 



E. Tagging Summary 

We summarize the tagging efficiencies we use in Ta- 
ble m For T and t decays, we only consider the fully 
hadronic decay channels. The efficiencies are used to 
determine the statistical uncertainties we give in our re- 
sults. For ratios of cross sections we assume that exper- 
imentalists have properly taken into account efficiencies 
to extract the appropriate cross section. 

TABLE I: Summary of the overall efficiencies used for esti- 
mating the number of events observed of the given fermion 
species. The overall efficiency is the tagging efficiency for the 
observed fermion decay mode times the BR to that final state. 



Channel 


Overall e/ 


Jet Rejection ej 


Z' — > /i^/i^ 


0.92 


n/a 


Z' T + T' 


0.08 


2.0 X 10"'' 


Z' bb 


0.36 


2.5 X 10"^ 


Z' tt 


0.075 


1.0 X 10""* 



F. Kinematic Cuts 

The values for the fake rates and efficiencies assume 
the ideal case of perfect alignment within the inner de- 
tector and no pile-up. Since tracking is only available 
within |r/| < 2.5, we do not include events in which one 
or both fermion tracks lie outside this region, and, unless 
otherwise stated, we require a minimum px > 20 GeV 
as needed for flavour tagging. Figure [T] (a) shows the 
invariant mass distribution in the bb final state for sev- 
eral models, the QCD bb backgrounds, and the light dijet 
backgrounds. 

For hadronic decays of r, 6 and t fermions, we found 
that a. pt > Q-ZMz' cut on the reconstructed momen- 
tum of the hadronic jets effectively reduced both the ir- 
reducible and dijet backgrounds as compared to the sig- 
nal. The events from the decay of a Z' tend towards 
a harder px distribution than the QCD backgrounds so 
that a larger percentage of Z' events pass this cut than 
the QCD backgrounds. Figure [T](b) shows the signal for 
several Z' models, QCD bb background, and the light di- 
jet background after applying this strong pT cut. This 
cut was also applied to muon final states when calculat- 
ing ratios of cross sections, for consistency. For deter- 
mination of the discovery limits we only apply a cut of 
Pt > 20 GeV to the final state muons. 

Fig. [T] (c) shows the reduction of the light dijet back- 
ground after taking into account detector efficiencies and 
fake rates. Figure [T] (c) shows that the application of 
appropriate flavour tagging algorithms reduces both the 
QCD bb and the light dijet backgrounds to a level where 
a meaningful measurement should be possible. 

As a final kinematic cut, we include only events within 
an appropriate invariant mass window around the reso- 
nance mass to improve the signal to background. Unless 
otherwise stated, we take this to be jM^/prime — Mj:j;\ < 
2.5r z/prime. This restricts the background events to 
the kinematic region directly under the resonance peak. 
However, this introduces an additional experimental un- 
certainty due to detector resolution, which smears out 
the resonance peak, effectively reducing the number of 
signal events in the peak. Because we use ratios of cross 
section measurements into bb, ti, fj,~^ ^~ , and r"'"r~ final 
states, which have different energy resolutions, the nu- 
merical value of the observables will shift. Thus, it will 
be important to understand detector resolution to accu- 
rately extract the underlying cross sections and coupling 
dependence. In addition, the reduction of the measured 
signal compared to background will increase the experi- 
mental errors. 

To properly account for detector resolution requires 
a realistic, detailed simulation for the specific particle 
identification algorithms being used. As discussed above, 
this subject is evolving rapidly. In addition, experimen- 
talists are constantly improving their understanding of 
the energy cahbration of the LHC detectors. To gauge 
the importance of detector resolution, we use estimates 
from recent detector studies and include them by ap- 
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FIG. 1: bb invariant mass distributions for QCD bb back- 
ground (solid, dark), light dijet background (solid, grey), and 
a Z' with Mz' = 2.0 TeV for the LRM (dot-dash), UUM 
(dot-dot-dash), and LH (dashed) models, (a) Only includes 
detector acceptance cuts of pr > 20 GeV and < 2.5. (b) 
Also includes pr > O.SMz' cut as described in the text but 
before including tagging efficiencies, (c) Includes the applica- 
tion of cuts, tagging efficiencies, and fake rejection rates given 
in Table ID 



plying Gaussian smearing to the final state momentum. 
Studies by the ATLAS collaboration expect between 3% 
and 5% energy resolution for TeV scale hadronic jets jll] 
(although the current ATLAS jet calibration gives 5% 
resolution [t^). This results in a resolution of ~3.5% for 
MijIj. This broadens the resonance, reducing the num- 
ber of measured signal events within the same invariant 
mass window. Signal significance may therefore be sig- 
nificantly reduced for narrow resonances. 

To see how detector resolution affects measurements, 
we include results for the ideal case of no smearing and 
for a more realistic case of 5% on b and t final states and 
3% on muon final states. For the bb case we found that for 
the chosen mass window, as expected, detector resolution 
has the greatest effect for the narrowest resonances (e.g. 
the Z'^) with virtually no effect on the broadest states 
(e.g. the ZIjum). 

The situation for the tt final state is less conclusive for 
the following reason. A recent ATLAS study [t^ gives a 
resolution of 4.6% for the tt invariant mass distribution 
of a narrow resonance. Using the same procedure used 
to gauge the importance of detector resolution for the 
bb channel we find that a usable signal can be measured 
in the ti channel for all models considered in this paper 
except for possibly Z,, and Z.^. This is consistent with 
the findings of Barger, Han and Walker [7lj. However, 
another ATLAS study 0| gives a resolution of 9-10% 
for a 1 TeV Z/p"'^'' with width T z/pr^^. /Mz/p.^,^. = 
3.3%. The range of these expectations demonstrates the 
difficulty in trying to predict the detector resolution for 
these measurements. 

Reconstruction of r+r^ final states are complicated 
due to the missing energy from the neutrinos in the t 
decays. LHC studies of t~^t~ final states with respect 
to Higgs searches have found a resolution on the recon- 
structed Higgs mass of Mh 10% [zl M,^- We do 
not incorporate this mass resolution in our results, but 
will refer to it when discussing our results for t~^t~ final 
states. 

Thus, while we do estimate the effects that detector 
energy resolution will have on the precision of our mea- 
surements, past experience shows that the experimental- 
ists eventually exceed initial expectations. In addition to 
this, the reduced signal to background caused by detec- 
tor resolution can be mitigated to some extent by higher 
luminosities and better identification efficiencies, which 
would improve the statistics. 

V. RESULTS 

A. Discovery Limits 

Leptonic final states offer the cleanest channel for the 
discovery of extra neutral gauge bosons due to the low 
backgrounds and clean identification [80l - [83 |. A very 
small number of dilepton events clustered in one or two 
bins of the invariant mass distribution would be taken 
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FIG. 2: Luminosity required for Z' discovery as a function 
of Mz' based on the observation of 5 /i^/i" events within 
the invariant mass window described in the text. The sets 
of curves are for an LHC centre of mass energy of -y/s = 
7 TeV and a/s = 14 TeV. ^From left to right, the models are 
ETC (dashed, dark), TC2 (solid, grey), UUM (dotted, dark), 
AFSLH (dashed, grey), SLH (solid, dark), SSM (dotted, grey) 
and LH (dot-dash, grey). 



as an obvious signal for new physics. To quantify this 
we consider two opposite sign leptons and impose kine- 
matic cuts of \r]i\ < 2.5 and pr, > 20 GeV to reflect 
detector acceptance. The criteria used was 5 events 
in the n'^ fJ-~ channel with a signal-over-background of 
at least 5 in an invariant mass window within ±1 bins 
of the resonance peak with the bin size as defined by 
AM = 24(0.625A/ + + 0.0056)^/2 Q^y ^^lere M is 
given in TeV j85|] . 

The integrated luminosity required to discover a Z' of 
a given mass in the dimuon channel is shown in Fig. [21 
We show curves for pp collisions with y/s = 7 TeV corre- 
sponding to the 2010-2011 LHC run and for ^/s = 14 TeV 
which corresponds to the LHC design energy. 

In Fig. |3]we show discovery limits for the various mod- 
els for several LHC benchmark energies and luminosities 
and compare them to discovery limits for the Fermilab pp 
collider. For the Tevatron we assume two cases; 1.3 fb^^ 
of integrated luminosity, as in Ref. [86|], and 8 fb~^ to 
estimate the reach for the full expected luminosity. We 
used similar detector acceptance and cuts as in Ref. [s^ : 
we impose a kinematic cut of px > 25 GeV and consider 
events within two regions of pseudorapidity - where both 
leptons satisfy \ri\ < 1.1, and where one lepton satisfies 
\r]\ < 1.1 and the other satisfies 1.2 < \r]\ < 2.0. As 
well, we consider only events within an invariant mass 
window of \Mz' - Mi+i- \ = ±10%Mz'. We use the dis- 
covery criteria of 5 observed dilepton events as with the 



1.96 TeV - 1.3 fb' 



1.96TeV-8.0fb' 



7 TeV -100 pb' 



7TeV- 1 fb 



14TeV- 1 fb" 



14TeV- lOfb 



14TeV- 100 fb 




FIG. 3: Discovery reach at benchmark luminosities for Teva- 
tron and LHC (both early and design) energies. 



LHC study. We note that this differs from that used by 
the CDF collaboration to obtain the current direct limits 
in Refs. 86] and Figure [3] shows that the 2010-2011 
LHC run will be able to roughly double the discovery 
reach of the Tevatron. 

In addition, we considered the possibility that a Z' 
with a preferential coupling to third generation fermions 
may be discovered first in the r+r" final state. For the r 
hadronic decay mode, none of the models we considered 
would be first visible in the r final state for a value of 
mixing angle, 0, that would give rise to a narrow reso- 
nance. At extreme values of the mixing angles in models 
without generation universality, the Z' begins to decou- 
ple from the first two generation fermions, reducing the 
qq — Z' coupling, and hence the production of the Z'. 
For these models, a Z' that would be first observed in 
the T+T^ final state would have a width larger than 10% 
of the Z' mass, resulting in a very broad peak, unlikely 
to be observed. However, as described by Holdom [t^I, if 
other T decay modes could be used, it may be possible to 
first observe a Z' in the t"*"t~ final state for some regions 
of parameter space. 

Constraints from electroweak precision data are more 
stringent than the discovery limits from the Tevatron for 
some models, as seen in Refs. [l^ H^. However, it 
is clear that the 7 TeV LHC run will improve the limits 
for models with universal couplings, and the full 14 TeV 
run should improve the limits for non-universal models 
for reasonable values of the mixing angle, (j). 
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B. Model Discrimination using tt and bb to /i fi 
Production Ratios 



The primary goal of this paper is to explore the use 
of third generation ferniions to distinguish between mod- 
els of extra neutral gauge bosons. We start with ratios 
of tt and bb to At'^M" cross sections, and expand on our 
previous study [7[ • We are particularly interested in mod- 
els with non-universal couplings - specifically, the UUM, 
ETC and TC2 models - as the R 



and Rb/^ ratios 



defined below produce results that are quite distinctive 
from the models considered in Ref. Q- 



Rt/n and Rb/p, are defined by: 



Rt/p. = 



a{pp ^Z' ^ tt) 3Kt {Li + 



6/m 



a{pp - 
a{pp' 



Z' 



bb) 



{LI 



3Kb{Ll+Rl) 



a{pp 



) 



(10) 



(11) 



where Lj and Rf are the left and right handed fermion 
couplings to the Z' and the K factors incorporate the 
QCD and QED NLO correction factors [H [Hj. The use 
of the Rt/fi and i?f,/^ ratios has the benefit of reducing 
the contributions from uncertainties in the parton dis- 
tribution functions as the initial state qq couplings and 
pdf's cancel in the ratio. 

To obtain our results, we assumed that a Z' has been 
discovered and its mass and width measured [13, US HO, 
[9l| such that the appropriate Mqq and pT cuts described 
above can be applied. We calculated the expected num- 
ber of events and statistical errors for signal plus back- 
ground for a given integrated luminosity and particle 
identification efficiencies, e^+^-, ^bin and eti from Ta- 
ble |T1 The expected number of SM QCD background 
events were subtracted from the total events to give the 
predicted number of signal events. We did not include 
systematic uncertainties arising from uncertainties in the 
luminosity and identification efficiencies. 

Our results for Rb/^ and Rt/^i are shown in Fig. S] (a) 
and (b), where Fig. |4] (b) expands the view to show 
the distinctive measurements that would be observed 
for models with preferential couplings. We show la 
statistical errors based on an integrated luminosity of 
L = 100 fb~^. The measurements for models with pref- 
erential couplings and mixing parameters that result in 
large couplings for third generation quarks are quite dis- 
tinctive, making this a possible hallmark measurement 
for these models. 

In Fig. [5] we show Rb/^ and i?t/^ including a detec- 
tor resolution of 5% for b and t final states, and 3% for 
muons. Comparing Fig. [S] with Fig. 31 it is clear that 
the effect of resolution is small for broad models such 
as ETC, TC2 and UUM. On the other hand, the effect 
is much more significant for the narrow models such as 
the E6, SLH, and AFSLH models, affecting both the 
values and expected statistical uncertainty. There are 
two approaches to account for resolution effects in these 
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FIG. 4: Measurements of the ratio of tt and bb to 
cross sections within 2.5r of Mz'- Black bars correspond to 
expected la statistical uncertainties for Mz' = 1.5 TeV and 
grey bars to Mz' = 2.5 TeV. (b) Expands the scales to include 
models with generation and family dependent couplings. 



measurements. The experimentalists can try to disentan- 
gle resolution effects from the underlying cross section to 
give the "theoretical" cross section or they can compare 
to Monte Carlo simulations that include detector reso- 
lution. In any case, the sensitivity of measurements for 
narrow models to detector resolution will lead to system- 
atic uncertainties in the measurements that need to be 
taken into account. 

For some models, including models with non-universal 
couplings, the Eq model, and the LR Symmetric mod- 
els, the predictions for Rb/p, and Rt/fj. are dependent on 
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FIG. 5: Measurements of the ratio of tt and bb to 
sections within 2.5r of Mz' including detector resolution of 
5% for b and t final states, and 3% for muons. The labelling 
is as in Fig. [l] 



tions for some models overlap in the Rt/fj. — Rb/^i space as 
shown in Fig. |6l Consequently, other measurements will 
be needed to distinguish between models for these values 
of the model parameters. 

It is important to note that the region of overlap in 
the Rt/ ^1 — Rb/fi plane between models with and without 
universal couplings occurs for parameter values that are 
not of particular interest. In the case of the UUM model, 
the overlap occurs for parameter values where leptons 
have preferential couplings and in the case of the ETC 
and TC2 models the overlap occurs when the first two 
generations of fermions have preferential couplings. Since 
these models are constructed such that the top quark 
plays a role in EWSB, one would not expect their mixing 
angles to take values in the overlap region. 

A final observation is that the UUM and ETC mod- 
els are indistinguishable using measurements of i?t/^ and 
Rb/n for any value of mixing parameter. However, in 
this case, the ratio of tau to muon events at the LHC 
will discriminate generation dependent couplings. This 
is examined in the next subsection. 



C. Extracting Mixing Parameters Using Rr/p 



the mixing angle between subgroups in the model. This 
is illustrated in Fig. |6] which plots the Rt/^ — Rb/fj, ra- 
tios while varying the mixing parameters. To obtain 
these curves, a 

Z/prune 

mass of 1.5 TeV was used, al- 
though the results are not very sensitive to M^/primc, 
and the following ranges for mixing angles were used: 
for the LR Symmetric model, both with standard and 
alternate isospin assignments, the mixing parameter is 
constrained by 0.55 < (gR/gL) < 1 123]; for the UUM, 
(j) is constrained by 0.22 < sm(j) < 0.99 |9^; and for 
the remaining models no specific limits could be found in 
the literature that were not directly tied to the mass of 
the Z' . Depending on the mixing parameter, the predic- 



Measuring generation universality will be an important 
step in distinguishing or ruling out TC2 and ETC type 
models. The simplest and cleanest way to measure the 
level of universality is to determine the ratio of the Z' 
decays to r-leptons and to muons or electrons which can 
be found by measuring <7{pp — > Z' — > t~^t~ ) / a{pp — >■ 
Z' fi+fi-). 

We also considered a measurement of the ratio of t- to 
c-quark cross sections but found that such a measurement 
does not appear promising due to low tagging efficiencies 
of the charm quark, and indistinguishability of the charm 
quark from the light jet backgrounds. 

We therefore restrict ourselves to measurement of the 
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FIG. 7: for models with universal and non-universal 

fermion couplings. The ratios are obtained for tan and muon 
events satisfying \Mj,f — Mz'\ < 2.5Tz', along with the other 
cuts described in the text. The dark bars are for Mz' ~ 
1.5 TeV and the grey bars for Mz' = 2.5 TeV (for sin<^TC2 = 
0.5 and sincjiETC ~ 0.9). 



Z' 



ratio i?r//j = cr(pp Z' ^ T^T^)/ij{pp 

shown in Fig. [T] To obtain these results we 
imposed a cut on the invariant mass of the final state 
fermions of \Mz' — -^//l < 2.5rz' in addition to impos- 
ing a requirement ofpT/ > O.SA/^' and |?7/| < 2.5. These 
results and those shown in Fig. [5] only include the sta- 
tistical errors for a semi-idealized detector with perfect 
energy resolution. We will comment on this below. 

It is clear from Fig. [7] that models with generation uni- 
versality will yield measurements of Rt/^ — 1 with rea- 
sonable precision. In contrast, models with generation 
dependent couplings show a large, measurable variation 
from unity. The dependence of Rr/^ on the mixing angle 
between the gauge groups of the theory are given by: 



ETC 



^iz'J,^c-^^^^^^-) 



cx tan"* 



'ETC 



and 



T{Z^ 



TC2 



r(z^ 



TC2 



^ OC COf* <PtC2 



(12) 



(13) 



Given the fundamental nature of the gauge group mix- 
ing angle, measuring its precise value would be vital in- 
put into constructing the Lagrangian of the underlying 
theory. In Fig. |5] we show how well such a measurement 
can be made for the TC2 and ETC models using Rr/^ 
assuming Mz' = 1.5 TeV and L = 100 fb~^ for the semi- 
idealized detector. In these plots the x-axis corresponds 
to the assumed value of the mixing parameter and the 
y-axis corresponds to the measured value of the mixing 
parameter with the spread in the vertical direction cor- 
responding to 1- and 2-(T limits of R-r/f^ for the input pa- 
rameter value and measured value. The parameter range 



corresponds to the range where the Z' width is less than 
10% of the Z' mass. As before, we include the back- 
grounds in estimating the statistical errors and impose 
the same kinematic cuts as before. These limits could be 
further constrained by including more observables, such 
as Rb/f_i and Rt/^ into the fit. 

Figure [7] indicates that, for typical parameter values, 
the non-universal models predict values for Rr/p. that 
are very distinct from the value of Rr/fj, — 1 expected 
for models with generation universality. As pointed out, 
these results do not include finite detector resolution and 
in addition, measurements of t^t~ final states have the 
additional complication of missing energy due to neutri- 
nos in the r decays. However, the non-universal ^/p''*™eig 
have large decay widths so that including the measure- 
ment resolution [tgI . [tsI [79j will only have a small effect 
on the results for non- universal z/p"™^'s. On the other 
hand, some of the generation universal Z^p""^" mod- 
els are relatively narrow so that the resonances will be 
smeared out in the r+r" final state and some care will 
have to be taken in choosing appropriate invariant mass 
windows. We note that both ATLAS ^ and CMS [H 
have recently measured (j{pp Z^) x i3(Z° — r+r~) 
with errors of roughly 20% and 10% respectively in the 
ThTi modes which bodes well for a measurement of Rr/^- 

We conclude that, at worse, a very crude measure- 
ment of Rt/^ could signal generation non-universality. 
The mixing parameter measurement capability shown in 
Fig. [5] is an idealization but, given the recent ATLAS 
and CMS measurements of Z° — > r+r^, we are opti- 
mistic that a measurement will be possible to constrain 
the relevant mixing angle. 



D. Forward-Backward Asymmetry with Heavy 
Quark Final States 

The ability to identify h- and t-quarks at the LHC of- 
fers the possibility of using forward-backward asymme- 
tries [Aps) in heavy quark final states to assist in the 
determination of individual fermion couplings to a Z' ^ . 
As will be seen below, ApB has a different dependence 
on the Z'-fermion couplings than the Z' production cross 
section and Z' width. 

Generally, a forward event is defined by the decay an- 
gle of the outgoing fermion relative to the direction of the 
interacting quark in the Drell-Yan annihilation. For pp 
collisions at the LHC, there is an ambiguity in determin- 
ing the direction of the quark, where it is impossible to 
tell on an event- by- event basis whether the Z' is boosted 
in the direction of the quark or anti-quark. Because the 
momentum distributions are harder for valence quarks 
than for sea anti-quarks, this ambiguity can be resolved 
to a certain extent by assuming that the Z' boost direc- 
tion is the same as the quark direction [96j . In the central 
region of Z' rapidity, the quark and anti-quark momenta 
are more evenly balanced, and the correct and incorrect 
assignments are canceled in an asymmetry measurement. 
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FIG. 8: LHC measurement capability of the mixing param- 
eter using R^i^ for (a) TC2 and (b) ETC. The horizontal 
axis corresponds to the input parameter value and the ver- 
ticle axis corresponds to the extracted value with 1- (grey) 
and 2- (black) a limits. Mz' ~ 1.5 TeV and an integrated 
luminosity of 100 fb~^ was assumed in these plots with the 
following kinematic cuts; j?7j < 2.5, \AIi+i- — Mz'\ < 2.hT z' , 
Pt < O.SMz'. 



In a recent paper [g] , we suggested a simple method of 
performing the forward-backward asymmetry measure- 
ment by using the direct pseudorapidity measurements 
of the final state particles. This method differs from 
the traditional definition given in 0, |93, ll^ • It can be 
shown that a "forward" event is one in which \rif \ > \rijr\ 



in the lab frame, and vice-versa for a "backward" event, 
assuming that the Z' is boosted by the quark as in the 
traditional definition. Using these forms for forward and 
backward events, the forward-backward asymmetry of 
the signal is given by: 



A 



FB — 



J m 


- B{y)dy 


J F{y) 


+ B{y)dy 
1 













(14) 



where F(y) is the number of forward events and B{y) is 
the number of backward events for a given Z' rapidity, 
y. The approximation is a representation of the on-peak 
contribution to the Aps that more clearly indicates the 
coupling dependencies. Lf and Rf are the left and right 
handed couplings of the Z' to the fermions, and are 
the integrated symmetric and anti-symmetric combina- 
tions of the parton distribution functions. 

This method for finding the ApB has the advantage 
of being very straightforward and clean. It simply relies 
on counting events with \rif \ > It^j] and \rif \ < \rij\. No 
calculation of the centre-of-mass scattering angle or Z' 
rapidity is required. In Fig. |9l we show the ApB distri- 
bution as a function of invariant mass of the bb and tt 
final states for several representative models with a Z' 
mass of 1.5 TeV. It should be noted that it is unlikely 
there will be sufficient statistics to make this measure- 
ment except on the resonance. 

As in all measurements involving third generation 
fermions, the challenge is extracting the events of interest 
from a large standard model background and accumu- 
lating sufficient statistics to make a meaningful measure- 
ment. The QCD backgrounds for these measurements are 
forward-backward symmetric at tree level, which should 
allow a heavy quark ApB to be sensitive to the presence 
of a Z'. Rather than subtracting out the backgrounds we 
include the totals of signal plus backgrounds in the for- 
ward and backward regions. Since the backgrounds are 
symmetric, they do not contribute to the numerator, but 
do for the denominator. This has the effect of making the 
magnitude of the asymmetry smaller but also has the ad- 
vantage of reducing the statistical errors. The statistical 
errors include contributions from both QCD heavy quark 
backgrounds and the reduced light dijet background. 

Figure [TO] shows the expected results for ApB in the 
bb and tt channels for several models assuming Mz' — 
1.5 TeV, L = 100 fb~^, and the same kinematic cuts de- 
scribed previously, not including detector resolution ef- 
fects. While large uncertainties are apparent for some 
models, a reasonable measurement can still be expected 
for most models, including the Left- Right Symmetric 
Model, the various Little Higgs models, and models with 
non-universal couplings. Table HIl gives the corresponding 
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FIG. 9: (a) Afg and (b) Apg of the Z' signal as a function of 
invariant mass for a Z' with AIz' = 1.5 TeV. Cuts employed 
include pt > O.SMz' and |r;| < 2.5. These figures do not 
include backgrounds, and only examine the asymmetry dis- 
tribution of pp — >■ j/Z^ /Z' — >■ qq. Models follow the same for- 
mat as Fig.[l ETC (dashed, dark), TC2 (solid, grey), UUM 
(dotted, dark), AFSLH (dashed, grey), SLH (solid, dark), 
SSM (dotted, grey) and LH (dot-dash, grey). The results are 
for the "ideal" case that doesn't take into account detector 
resolution. 



numerical values for A^g'' and the statistical uncertain- 
ties that can be expected for 100 fb~^ integrated lumi- 
nosity. 

Fig. llll shows similar results but this time including the 
detector resolution as described earlier. In this figure we 
assumed an integrated luminosity of 300 fb~^ to improve 
the statistics. Taking into account detector resolution 
both reduces the statistics and shifts the observed val- 
ues for ApB- While ApB still has some resolving power 
with bb final states it is not clear how useful the tt final 
states will be, due to low statistics. It is possible that if 
some effort were made to disentangle detector resolution 
from the underlying cross section these results could be 
improved. Nevertheless, we expect that ApB measure- 
ments would help constrain the Z'-fermion couplings as 
part of a global fit. 




FIG. 10: A^l^B (black bars) and A^J^g (grey bars) for a Z' with 
a mass of 1.5 TeV. Statistical errors include contributions 
from QCD backgrounds and light dijets assuming 100 fb^^ 
luminosity. 
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FIG. 11: A''pg (black bars) and Apg (grey bars) including 
detector resolution as described in the text. Statistical errors 
include contributions from QCD backgrounds and light dijets 
assuming 300 fb~^ luminosity. 



VI. SUMMARY 

Extra neutral gauge bosons are a hallmark of many 
models of physics beyond the standard model and may 
be discovered early in the LHC program. In this pa- 
per, we expanded upon our previous phenomenological 
study by focusing on measurements using third genera- 
tion fermions and exploring models with non-universal 
fermion-Z' couplings. We first gave an update of Z' 
discovery limits including non-universal coupling mod- 
els and the LHC energy and planned luminosity for the 
2010-2011 run. We found that this run will roughly dou- 
ble the mass reach of the Fermilab Tevatron. 

The main focus of this paper was to explore the useful- 
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TABLE II: Afb values for b- and i-quark final states with corresponding statistical uncertainties, assuming Mz' ~ 1.5 TeV 
and L = 100 fb"\ Cuts include pr > 0.3Mz> GeV, \r]t,b\ < 2.5, within \ AMgg - Mz' \ < 2.bVz'. The first set of data, labelled 
as "ideal detector", do not include detector resolution smearing for the final state fermions, while the second set of data does 
include the effects of detector resolution. 



Model 


4* 

^FB 


± 


5A'fb 




SA^FB 


A^FB i 


5A'fb 




SA'-FB 








ideal detector 






including 


detector resolution 




Ed X 


0.00 


± 


0.02 


0.060 ± 


0.011 


0.00 ± 


0.03 


0.017 ± 


0.016 


Ee tp 


0.00 


± 


0.03 


0.000 ± 


0.016 


0.00 ± 


0.05 


0.00 ± 


0.02 


Ee n 


0.00 


± 


0.03 


0.013 ± 


0.014 


0.00 ± 


0.04 


0.00 ± 


0.02 


LRM 


0.097 


± 


0.016 


0.181 ± 


0.007 


0.04 ± 


0.02 


0.075 ± 


0.011 


ALRM 


0.146 


± 


0.019 


0.032 ± 


0.013 


0.04 ± 


0.03 


0.006 ± 


0.016 


UUM 


0.144 


± 


0.007 


0.189 ± 


0.004 


0.073 ± 


0.011 


0.094 ± 


0.006 


SSM 


0.049 


± 


0.014 


0.112 ± 


0.007 


0.02 ± 


0.02 


0.044 ± 


0.010 


TC2 


0.069 


± 


0.017 


0.027 ± 


0.010 


0.02 ± 


0.02 


0.01 ± 


0.01 


LH 


0.114 


± 


0.012 


0.158 ± 


0.006 


0.047 ± 


0.018 


0.063 ± 


0.009 


SLH 


0.086 


± 


0.018 


0.136 ± 


0.009 


0.03 ± 


0.03 


0.042 ± 


0.014 


AFSLH 


0.079 


± 


0.019 


0.124 ± 


0.010 


0.02 ± 


0.03 


0.035 ± 


0.015 


3-3-1 (2U ID) 


0.050 


± 


0.018 


0.083 ± 


0.010 


0.02 ± 


0.03 


0.023 ± 


0.014 


ETC 


0.058 


± 


0.011 


0.086 ± 


0.006 


0.022 ± 


0.015 


0.031 ± 


0.008 



ness of third generation fermions in studying extra neu- 
tral gauge bosons. We found that it should be possible 
to measure the decays of a moderately heavy Z' to third 
generation quarks, depending on the capability of the 
experiments to reject against the light jet background. 
Such measurements would prove to be very effective at 
distinguishing between different models of Z' . The mea- 
surement of the ratio of r to /i decays of the Z' should 
be very effective in testing generation universality of the 
Z' . We also studied using the forward backward asym- 
metry of h- and i-quarks from Z' decays to distinguish 
models. For a Z' of moderate mass these measurements 
could help distinguish between different models. More 
importantly, they can contribute valuable input for mea- 
suring Z'-fermion couplings, and should be an integral 
piece of a global fit with this goal. 

An important source of uncertainty in these measure- 
ments is the degradation of the heavy fermion final state 
signals due to detector resolution. It appears to be man- 
ageable for hh and ti final states for broad reso- 
nances such as in the TC2 and ETC models, but becomes 
increasingly important for narrow resonances such as the 



E6 models. Now that ATLAS and CMS have demon- 
strated the ability to measure (j{pp — >■ Z^) x B{Z^ — >• 
T"'"r~) we are optimistic that R-r/fi can be measured well 
enough to distinguish between typical generation non- 
universal models and generation universal models. Fi- 
nally, we found that detector resolution effects can de- 
grade ApB measurements quite significantly, in particu- 
lar for ti final states. We believe that measurements us- 
ing 3rd generation fermions look promising but clearly, 
a more careful detector level study is needed to properly 
take into account detector resolution and other compli- 
cations of a real detector. 
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